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, sources have allowed the determination of the rest-frame 2 — 8 keV luminosity 



Recent deep X-ray surveys combined with spectroscopic identification of the 



as a function of redshift. In addition, an analysis of the HEAOl A2 2 — 10 keV 



full-sky map of the X-ray background (XRB) reveals clustering on the scale of 
several degrees. Combining these two results in the context of the currently 
favored AC DM cosmological model implies an average X-ray bias factor, b x , of 
b\ = 1.12 ± 0.33, i.e., b x = 1.06 ± 0.16. These error estimates include only 
statistical error; the systematic error sources, while comparable, appear to be 
sub-dominant. This result is in contrast to the large biases of some previous 
estimates and is more in line with current estimates of the optical bias of 
^ ! galaxies. 



Subject headings: large-scale structure of the universe — X-rays: galaxies — X- 
rays: general 



1. Introduction 

An important test of any cosmological model is that it be consistent with the observed 
distribution of matter in the universe. Since our primary knowledge of this distribution 
comes from observations of galaxies, it is essential to understand the extent to which galaxies 
trace the matter density. This relationship is usually quantified by a bias factor which 
relates fluctuations in the galaxies to those in the dark matter. It is complicated by the 
fact that the relation between the luminosity of a galaxy or groups of galaxies and the 
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underlying distribution of matter can depend on the type of galaxy, the spectral band of 
the observation, the redshift z, and the scale length on which the comparison is made. 
However, such complications are also opportunities in that models of galaxy formation must 
successfully reproduce these differences. 

The standard definition of the bias factor, 6, is the ratio of the fractional galaxy density 
fluctuations to the fractional matter density fluctuations, i.e., 

b = *hlpL (i-i) 

bp/p 

where p g is the mean density of galaxies, p is the mean density of matter and b indicates 
the rms fluctuations of the densities about these means. If galaxies formed early (z > 1), 
as appears to be the case (e.g., Ellis 1997), then there are good reasons to expect that, for 
linear density perturbations (i.e., Sp/p <C 1) on large scales in the nearby [z < 1) universe, 
galaxies should be relatively unbiased (6 — > 1) tracers of the density field (Fry 1996, Tegmark 
& Peebles 1998); however, this assertion must be tested. 

Here we focus on determining the bias of the hard X-ray background (XRB), which 
is known to be dominated by distant [z < 2) active extragalactic galaxies and so provides 
a probe of the bias on large scales. The observed clustering of the XRB, when combined 
with what is known about the level of perturbations and the cosmological model from CMB 
observations, allows us to place relatively strong constraints on the X-ray bias. 

Previous determinations of X-ray bias have resulted in a wide range of values, 1 < b x < 7 
(see Barcons et al. 2000 and references therein). Some spread in the estimates is to be 
expected; e.g., at lower energies X-ray emission is dominated by clusters of galaxies, and 
so are expected to be as highly biased as clusters themselves (Bahcall & Soneira 1983). 
However, another major contribution to the uncertainty in the bias estimates is the lack 
of accurate determinations of the clustering of various X-ray sources. Two of the lower 
estimates of X-ray bias are from Treyer et al. (1998) who found that b x ~ 1 to 2 from a low 
order multipole analysis of the HEAOl A2 data set and Carrera et al. (1998) who found 
that the ratio of the X-ray bias of AGN to that of IRAS galaxies to be 0.8 < b x /bj < 1.7 
from ROSAT observations. The remaining uncertainty in these determinations arose from 
uncertainties in both the X-ray luminosity function (LF) and in the cosmological model. 
Knowledge of both of these has improved dramatically in the last year and this is largely 
responsible for the improved accuracy of the estimate of b x in this paper. 
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2. Clustering in the HEAOl A2 2 - 10 keV X-ray Map 

We recently presented evidence of large-scale clustering in the HEAOl A2 2—10 keV full- 
sky map of the hard XRB on angular scales of < 10° (Boughn, Crittenden, & Koehrsen 2002). 
Before computing the correlations, local sources of the X-ray background were removed from 
the map. The map was masked so as to eliminate strong, nearby X-ray sources with fluxes 
exceeding 3 x 10~ n erg s ^cmT 2 . In addition, all regions within 20° of the Galactic plane 
or within 30° of the Galactic center were masked. The map was also corrected for a linear 
time drift of the detectors, high Galactic latitude diffuse emission, emission from the local 
supercluster, and the Compton-Getting dipole. The latter components were fit to and then 
removed from the map. Without these cuts and corrections, the correlations are dominated 
by a few strong point sources and large-scale diffuse structures in the map. 
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Fig. 1. — The auto-correlation function of the HEAOl A2 map with bright sources and the 
Galactic plane removed. The dashed curve is the form expected from beam smearing due 
to the PSF of the map while the solid curve includes a contribution due to clustering in the 
XRB. 
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Figure 1 is a plot of the intensity angular correlation function (ACF) given by 

ACF(6) = J> - /)(/, - I) (2-1) 

where the sum is over all pairs of map pixels, i, j, separated by an angle 6, l{ is the intensity 
of the ith pixel, / is the mean intensity, and Ng is the number of pairs of pixels separated 
by 9. Photon shot noise only appears in the 8 = 0° bin and has been removed. The 
highly correlated error bars were determined from 1000 Monte Carlo trials in which the 
pixel intensity distribution was assumed to be Gaussian with the same ACF as in the figure. 
The dashed curve represents the expected functional form of the contribution to the ACF 
due to telescope beam smearing of a random distribution of uncorrelated sources normalized 
to the ACF(O) point. It represents the profile that is expected if there were no intrinsic 
correlations in the XRB. The point spread function (PSF) of the map is due to pixelization 
and to the finite telescope beam and was accurately determined from the profiles of 60 bright, 
nearby point sources. It is clear from Figure 1 that the XRB possesses intrinsic (i.e., not 
due to beam smearing) correlated structure out to angular scales of ~ 10°. Full details of 
the analysis are discussed in Boughn, Crittenden, & Koehrsen (2002). 



3. The 2 — 10 keV X-ray Luminosity Function 

In order to determine the X-ray bias factor b x from the measured ACF, it is essential to 
know from which redshifts the X-ray fluctuations originate; the underlying density fluctua- 
tions grow quickly, so it is important that they be compared to the X-ray fluctuations at the 
same redshifts. This requires understanding the contribution to the 2 — 10 keV X-ray LF as a 
function of redshift. However, the HEAOl A2 observations are total intensity measurements 
of the hard XRB with no information as to the fluxes or redshifts of individual sources, so 
we must infer the LF by other means. Recently the Chandra satellite has made possible 
large, faint hard X-ray surveys with measured redshifts. Cowie et al. (2003) and Steffen et 
al. (2003) have combined Chandra sources with brighter sources from ASCA (Akiyama et 
al. 2000) and ROSAT (Lehmann et al. 2001) to determine the redshift evolution of the 
2 — 8 keV LF with few assumptions about the character of the sources. The incompleteness 
uncertainty in the redshift dependence of the volume X-ray emissivity is estimated to be 
a factor < 2 at any redshift. The spectroscopically identified sources comprise 75% of the 
total 2 — 8 keV X-ray intensity; the dominant uncertainties result from the unknown redshift 
distribution of the unidentified sources. 
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Fig. 2. — The volume emissivity as a function of redshift. The local value (z = 0) is the 
measurement of Miyaji et al. (1994). The low redshift points are derived from the data of 
Steffen et al. (2003) while the high redshift [z > 1) come from Cowie et al. (2003). 
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Emissivity as a function of redshift, X x {z), is plotted in Figure 2. The 2—8 keV emissivity 
in the redshift range 1 < z < 4 is taken to be that estimated by Cowie et al. (2003) using 
ROSAT data. In the range 0.1 < z < 1.0, we use the emissivity implied from the luminosity 
function of Steffen et al. (2003). Finally, for z = 0, we use the value of the local emissivity 
from Miyaji et al. (1994). The models discussed below are based on a polynomial fit which 
passes through the data points; however, the results are largely independent of the details 
of the fitted function. 
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Observed 


Best 


Low z 


High z 


Ueda et al. 


i{z < 1) 


54%/58% 


57% 


67 % 


48 % 


52 % 


b x 




1.06 


0.86 


1.36 


1.12 



Table 1: Properties of four models of X-ray emissivity: fraction of the intensity arising from 
z < 1, f(z < 1)) and implied bias, b x . See text for details. 
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The HEAO data is band limited and the X-rays detected at high redshifts have larger 
rest frame energies than they do locally. Therefore, in order to apply K-corrections to 
the observed intensities and to transform from 2 — 8 keV emissivities to the 2 — 10 keV 
values appropriate for the HEAO map, we must make an assumption about the frequency 
and redshift dependence of the XRB. If the ISM column density in front of an AGN is 
large enough (Nh > 10 21 cm 2 ) the observed spectrum will be hardened by photo-electric 
absortion. At high redshifts, the rest frame energies of the detected X-rays are relatively 
large and the effect of photo-electric absorption is less. Therefore, for a given column density, 
sources at high redshift will appear softer than their low redshift counterparts. As a crude 
approximation of this affect we assumed a photon spectral index T(z) = 1.2 + 0.2z where 
dN/dE oc E~ r is the number spectrum of the photons of energy E. This roughly describes 
the redshift dependence of sources with an intrinsic spectral index of Y — 1.8 subject to 
photoelectric absorption by column densities of Nh ~ 10 22 — 10 23 cm" 2 . Furthermore, the 
flux weighted average spectral index is V — 1.4 as is observed. The difference between this 
model and one that assumes a constant spectral index of T(z) = 1.4 is not large in the 
sense that the intensity distributions, dl/dz, of the two models fall well within the range 
of the extreme models of Figure 3. It is also possible that higher redshift AGN are more 
heavily absorbed (e.g. Worsley et al. 2004); however, such effects are not included in some 
models of the XRB (e.g., Ueda et al. 2003). If the absorption levels were higher at higher 
redshifts, this would reduce the evolution of the spectral index, i.e., lessen it's dependence 
on redshift. The fact that we observe more energetic rest frame photons is balanced by the 
fact that the sources are more absorbed. While neither of the models (constant T(z) = 1.4 
or r(z) = 1.2 + 0.2,2) is likely to accurately describe the actual spectrum, the fact that the 
biases of these models are within a few percent of each other indicates that uncertainty in 
the redshift dependence of the spectral index of the XRB is not an important source of 
systematic error. 
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Fig. 3. — Four different models for the contribution to the x-ray luminosity as a function of 
redshift. The solid line is our best estimate given the volume emissivity. The long (short) 
dashed model results from pushing the emissivity distribution to lower (higher) redshifts. 
Finally, the dotted line represents the model of Ueda et al. (2003), which is not constrained 
below z — 0.1. 
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lt is straightforward to compute the intensity distribution, dl/dz, from X x (z) in the 
context of the ACDM cosmological model used by Cowie et al (2003) (fi m = 1/3, Q\ = 2/3, 
and H = 65 km s~ l Mpc~ 1 ). While this model is somewhat different from that currently 
favored by the WMAP satellite data (Spergel et al. 2003), dl jdz is a directly observable 
quantity that is independent of the cosmological model. The dl jdz resulting from our canon- 
ical emissivity model is given by the middle solid curve in Figure 3 where the normalization 
is fixed by the LF of Cowie et al. (2003) and Steffen et al. (2003). 

This profile implies that the bulk of the XRB arises at much lower redshifts than previ- 
ously thought (e.g., Comastri et al. 1996) as was first pointed out by Barger et al. (2001). 
For this particular model, 57% of the 2 — 10 keV background arises from redshifts less than 
1. This is in agreement with the recent observations of Barger et al. (2003) that indicate 
54% of the spectroscopically identified 2 — 8 keV intensity arises at z < 1. This increases to 
58% when photometric redshifts are included (Barger et al. 2002; Barger et al. 2003). The 
total integrated intensity of our canonical model, 5.6 x 10~ s erg s^cm^sr' 1 , lies between 
and is consistent with both the HE AO estimate of 5.3 x 10~ 8 erg s^cm^sr^ 1 (Marshall et 
al. 1980; Gruber et al. 1999) and that estimated from ASCA satellite data (Gendreau et 
al. 1995; Kushino et al. 2002), 6.4 x 10~ 8 erg s^cm^sr^ 1 . In any case, the current analysis 
only requires the functional form of dl/dz and not the overall normalization. Finally, the 
intensity weighted spectral index of the model, f = 1.40, is the same as that observed for 
the hard XRB (Marshall et al. 1980; Gendreau et al. 1995). 

In order to test the sensitivity of the implied X-ray bias to the LF, we consider three 
alternative models of dl/dz. The upper dashed curve in Figure 3 has been weighted to low 
z by squeezing (in redshift) the canonical emissivity by a factor of 0.8 while fixing the local 
emissivity to be the 1 a upper limit of Miyaji et al. (1994). This model is fairly extreme, 
as it overestimates the intensity coming from low redshifts. (See Table 1 for a summary 
of the properties of the various X-ray models). The lower dashed curve in Figure 3 was 
weighted to high z by stretching the canonical emissivity by a factor of 1.3 while fixing the 
local emissivity to be the la lower limit of Miyaji et al. (1994). This model significantly 
underestimates the intensity coming from z < 1. Finally, the dotted curve in Figure 3 is from 
the recent AGN synthesis model of Ueda et al. (2003). Unfortunately, their model of dl/dz 
did not extend below z — 0.1 and our results below depend somewhat on the behaviour 
assumed for low redshifts. 
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4. Matter Fluctuations in a ACDM Universe and X-ray Bias 

Given the X-ray luminosity function, the linear bias factor can be inferred from the 
cosmological model, but only if the time dependence and scale dependence of the bias are 
known. In our analysis we will assume both redshift and scale independence of the X-ray bias. 
Given our nominal intensity distribution, and assuming the standard ACDM cosmology, the 
dominant contribution to the ACF on angular scales of a few degrees comes from structures 
with redshifts, 0.03 < z < 0.5, which correspond to linear scales of from approximately 
10 Mpc to 200 Mpc. This is a strong indication that we are observing clustering in the 
linear regime and so can use the straightforward analysis of the growth of linear structures 
in a ACDM universe. 

Using the current WMAP ACDM parameters (Spergel et al. 2003), i.e., Q m = 0.27, 
Q\ = 0.73, and Hq = 71 km s~ l Mpc~~ 1 , it is straightforward to compute a projected matter 
ACF with the same redshift distribution as for the canonical model (e.g., Boughn, Crittenden, 
& Turok 1998). If our assumptions about the bias are correct, the intrinsic X-ray ACF should 
have the same shape as the normalized matter ACF, with a relative amplitude given by the 
square of b x , the X-ray bias factor. 

The observed ACF also contains components due to beam smearing of uncorrelated 
X-ray sources and photon shot noise, the latter of which is uncorrelated and, therefore, only 
contributes to the ACF at 9 — 0. Therefore, any fit to the full data set must include these 
three components. At 9 = 0° the ACF is dominated by beam smearing and photon shot 
noise while above 9 = 12° the signal to noise is small. The solid curve in Figure 1 is the two 
parameter, maximum likelihood fit to the data in the range 2.5° < 9 < 12°. The implied 
X-ray bias is b 2 = 1. 12 ± 0.33 (1 a error) or b x = 1.06±0.16 with a x 2 of 4.6 for 6 degrees of 
freedom. Since the distribution of errors in the ACF is to a good approximation Gaussian, 
the statistical error attached to b x 2 as well as the x 2 of the fit have the usual interpretations. 
The error indicated for b x represents the 68% confidence interval; however, this error is not 
Gaussian. The signal to noise of the data point at 5.2° is 4 a and a variety of fits (see below) 
of b x 2 to the ACF indicate statistical significances between 3 and 4 a. 

We performed a variety of other fits to the data to check the robustness of our estimate of 
b 2 . A three parameter fit to the data in the full interval (0° < 9 < 12°) gives b x = 0.96±0.16. 
A one parameter fit for the large angle correlations (5.2° < 9 < 12°), where the beam 
smearing component is nearly negligible, gives b x = 1.20 ± 0.14; even a fit to the single 
datum at 5.2° yields a consistent value of b x = 1.25 ± 0.16, though it is probably mildly 
contaminated by the beam smearing component. Following our previous work (Boughn, 
Crittenden, & Koehrsen 2002), we also modeled the clustering term as a power law, oc l/9 a , 
with 0.8 < a < 1.6. These fits varied in amplitude; however, at 9 — 4.5° all of the fits agreed 
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to within a few percent. Normalizing the model clustering ACF to this level implies a bias 
of b x = 1.06 ± 0.17, also consistent with our canonical fit. The reduced xV s f° r these fits are 
all ~ 1 and the fits are all consistent with each other. 

The process of fitting for large-scale, diffuse components and then removing them from 
the RE AO map, results in some attenuation of the ACF on angular scales > 10°. These 
factors were determined from the same Monte Carlo trials that were used to determine the 
statistical errors and the fits were adjusted accordingly. Even if these factors are ignored, 
the fit value of b x changes by only 3%. 

To evaluate the level of uncertainty due a systematic error in the intensity distribution of 
the XRB, the two "extreme" models of Figure 3 were also fit to the data in the 2.5° < 9 < 12° 
interval. The biases resulting from these two fits are b x = 0.85 ± 0.13 for the low z model 
and b x = 1.36 ±0.21 for the high z model. Since these models are somewhat exaggerated, we 
conclude that they represent lower and upper limits of systematic errors due to uncertainty 
in dl/dz. A fit to the Ueda et al. (2003) model indicated in Figure 3 results in a similar 
value of the bias, though the precise results depend on how the model is extended to low 
redshifts (0 < z < 0.1). If this model is extended so that the low z behavior is not allowed 
to fall below that implied by Miyaji et al. (1994), dl/dz = 2.7 x I0~ 8 erg s^cm^sr' 1 , 
then the fit value of b x becomes 1.12 ± 0.17, which is consistent with that implied by our 
canonical model. If instead, we use a linear extrapolation to low redshifts, the bias can be 
somewhat (~ 15%) higher, but the local emissivity of this model would be nearly 2a below 
that implied by Miyaji et al. (1994). 

The ACF on large angular scales is quite sensitive to the contribution of low z sources 
(roughly half the ACF at 9 = 4.5° is due to sources with z < 0.1), so any error in estimating 
the low redshift cutoff in dl/dz could affect the results dramatically. By masking sources 
stronger than 3 x 10~ n erg s _1 cm~ 2 we effectively truncate the intensity distribution at low 
redshifts. The truncated profiles were determined from the flux cut and the local luminosity 
function of Steffen et al. (2003). If the value of the flux cut is in error due to, for example, 
a difference in normalizations of the source catalog used to make the cuts (Piccinotti et al. 
1982) and the Steffen et al. luminosity function, then this would affect the cutoff redshift and 
would be translated to an error in the predicted ACF. In the extreme limit of no flux cut, i.e., 
no truncation of the dl/dz profile, the implied X-ray bias is b x = 0.90. In the other extreme, 
i.e., a flux cutoff of 1 x 10~ n erg s^ 1 cm -2 , the implied bias is b x = 1.13. Therefore, it is 
unlikely that inaccuracy in characterizing the flux cut is the source of significant systematic 
error. 

Potentially more problematic is the redshift distribution of the unresolved component 
of the XRB. Worsley et al. (2004) found that above 7 keV only ~ 50% of the XRB is 



-14- 



resolved; although, this conclusion must be tempered somewhat by the fact that the brightest 
sources they considered (in the Lockman Hole) have fluxes of ~ 10~ 13 erg s~ l cm' 2 . Sources 
brighter than this contribute to the whole-sky XRB and they conclude that the true resolved 
fraction may be 10 to 20% higher. Even though only ~ 20% of the counts in the HEAO 
passband comes from photons with energies above 7 keV, an unresolved component can still 
significantly affect the estimate of the bias. As a pessimistic case, we ignore the bright source 
correction and assume a 30% unresolved component below 5 keV. In this case roughly 1/3 
of the 2-10 keV XRB is unresolved. If this unresolved component is distributed in redshift 
like the resolved component, then there is no change in the implied bias. On the other hand, 
if the unresolved component is entirely due to sources at high redshift where it does not 
contribute to the ACF signal, then the implied bias will be 50% higher than our canonical 
value. If instead the unresolved component is due to sources at low redshift, z < 1, then the 
implied bias will be 20% lower than our canonical value. These fall somewhat outside our 
two "extreme" values in Table 1 and so provide a caveat to those estimates of the limits of 
systematic error. However, if only half of the unresolved component is located at high (low) 
redshifts and the other half is distributed like the resolved component, then the implied bias 
is only 20% (11%) higher (lower) than our canonical value, well within the limits of Table 1. 

It is difficult to quantify all possible systematic errors; however, considering that the 
above "extremes" result in errors of the same order as the statistical error in the fit, we 
conclude that the total systematic error is no larger than the statistical error quoted. 

5. Discussion 

We have determined the X-ray bias of the hard XRB assuming it is time (i.e., redshift) 
and scale independent. These assumptions are probably quite reasonable since the mean 
redshift weighting of the X-ray ACF is quite low, z ~ 0.1, and the linear scales probed by 
the ACF are quite large (10 Mpc to 200 Mpc). Even if these assumptions are violated to 
some extent, b x can still be interpreted as an 'average' X-ray bias. There are several types 
of sources that contribute to the XRB, including quasars, Seyfert galaxies, LINERS, and 
clusters of galaxies, and the implied value of the bias must be considered to be an average 
over all these sources. However, the dominant contribution to the XRB is most likely to be 
moderately active AGN (Cowie et al. 2003), so b x should be representative of the bulk of 
the sources of the XRB. 

With these caveats in mind, we find an X-ray bias of b x = 1.12±0.33, i.e., b x = 1.06±0.16 
(statistical error only). This error includes photon shot noise, fluctuations in the XRB 
from beam smearing, and the clustering of the XRB itself. The fits of b x for two extreme 
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models of dl/dz indicate that the uncertainty due to our ignorance of the X-ray luminosity 
function is likely less than the statistical error. Other possible sources of systematic error 
also seem small. We conclude that the hard XRB is a largely unbiased tracer of the matter 
distribution on large scales. This is consistent with current models of large-scale, late time 
galaxy biasing (Benson et al. 2000; Tegmark & Peebles 1998). In addition, the latest studies 
of the clustering of ~ galaxies on ~ 100 Mpc scales indicates that these objects are also 
unbiased tracers of matter. Verde et al. (2002) found that, on scales of ~ 7 to ~ 40 Mpc, 
b = 1.04±0.11 for 1.9 galaxies in the 2dF survey with a mean redshift of z = 0.17. Using 
a different analysis of the same data, Lahav et al. (2002) found that b = 1.20 ± 0.11 on 
scales of ~ 20 to ~ 150 Mpc. Both of these results are consistent with early findings from 
the SDSS and 2MASS surveys that imply linear bias factors on the order of unity (Tegmark 
et al. 2002; Miller et al. 2003). It should not be surprising that the XRB and galaxy biases 
are similar since L* galaxies are closely associated with the moderately active AGN that 
comprise the bulk of the hard XRB (e.g., Barger et al. 2003; Miller et al. 2003). 

If these estimates are accurate, then the X-ray bias factor in the linear regime is now 
much better determined. The hard XRB background appears to be an excellent tracer of the 
large-scale distribution of matter, making it a useful tool for understanding the evolution 
of structure in the universe. One example of the importance of determining galaxy biases 
(and indeed the driving motivation for this work) is to aid in the interpretation of recent 
detections of correlations of galaxies with the cosmic microwave background (CMB). We 
(Boughn & Crittenden 2004) have detected a correlation of the 2 — 10 keV XRB with 
WMAP satellite map of the cosmic microwave background (Bennett et al. 2003), and there 
have been correlations observed with a number of other galaxy surveys (Nolta et al. 2003; 
Scranton et al. 2003; Fosalba, Gaztanaga,& Castander 2003; Afshordi, Loh, & Strauss 
2003). These correlations have been interpreted as the detection of the integrated Sachs- 
Wolfe (ISW) effect (Sachs & Wolfe 1967). If confirmed, they would constitute an important 
test of the AC DM cosmological model and provide further evidence of the existence of a 
substantial amount of "dark energy" in the universe (Crittenden & Turok 1996). 

We would like to acknowledge Keith Jahoda who is responsible for constructing the 
HEAOl A2 X-ray map and who provided us with several data-handling programs. We also 
thank Greg Koehrsen for noise analysis programs. RC acknowledges financial support from 
a PPARC AF fellowship. 
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